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Abstract. We have developed a detailed standard chemi- 
cal evolution model to study the evolution of all the chem- 
ical elements up to the iron peak in the solar vicinity. We 
consider that the Galaxy was formed through two episodes 
of exponentially decreasing infall, out of extragalactic gas. 
In a first infall episode, with a duration of ~ 1 Gyr, the 
halo and the thick disk were assembled out of primordial 
gas, while the thin disk formed in a second episode of in- 
fall of slightly enriched extragalactic gas, with much longer 
timescale. The model nicely reproduces the main observa- 
tional constraints of the solar neighborhood, and the cal- 
culated elemental abundances at the time of the solar birth 
are in excellent agreement with the solar abundances. By 
the inclusion of metallicity dependent yields for the whole 
range of stellar masses we follow the evolution of 76 iso- 
topes of all the chemical elements between hydrogen and 
zinc. Those results are confronted with a large and recent 
body of observational data, and we discuss in detail the 
implications for stellar nucleosynthesis. 

Key words: Nuclear reactions, nucleosynthesis, abun- 
dances - Stars: abundances - Galaxy: abundances - 
Galaxy: evolution - Galaxy: general 



1. Introduction 

The chemical evolution of our Galaxy has been exten- 
sively studied during the last years. Along the last decade 
a great deal of new data have become available, reveal- 
ing in principle the chemical history of the Milky Way 
and building up a set of observational constraints that 
have to be fulfilled by any successful theoretical model. 
These observational results concern in partic ular the age- 
metallicity relation (Meussingeret al. 1991 : Edvardsson 
et al. |l 993t Rocha-Pinto et al. |2000 ), star metallicity dis- 



tribution (Wyse fc Gi lmore 1995 ; Rocha-Pinto & Maciel 
1996 ; J0rgensen |2000 ), abundance ratios of an increas- 



ing number of chemical elements, both in halo a nd lo- 
cal disk stars (between many others, Sneden et al. 1994 
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McWilliam et al. |1995| ; Ryan et al. [L99§ Chen et al. pOOOj ), 
and radial abundance gradients (see Henry & Worthcy 



1999 for a review). 



These observational efforts have been accompanied by 
the publication of several theoretical models which try to 
interpret the data. Most of that work are standard open 
chemical evolution models (Carigi 1994; Giovanoli & Tosi 
1995| ; Prantzos & Aubert |1995| ; Timmes et al |1995| (here- 
aft er TW W1995); Chiappini et al. |l997j |1999| ; Tho mas et 
al. |1998| ; Chang et al. |l999| ; Goswami & Prantzos |2000 
(hereafter GP2000), etc.) in which the Galaxy is assem- 
bled by infall of extragalactic gas and that use different 
prescriptions for the main ingredients, i.e. infall rate and 
composition of the accreted material, initial mass func- 
tion (IMF), star formation rate (SFR) and stellar yields. 
Chemodynamical models have also been published (Stein 
Samland et al 



metz & Mller 1994 



1997; Samland 1998 



Berczik 1999), but due to the numerical complexity they 
have to introduce several approximations in the treatment 
of stellar lifetimes and yields. 



As pointed out by GP2000, while the prescriptions for 
the IMF, the SFR and the infall rate are just empirical 
recipes due to our poor understanding of the physical pro- 
cesses involved, the stellar yields can be determined from 
first principles, as they rely upon the much better known 
theory of stellar evolution. However, there still remain se- 
vere uncertainties in the stellar yields, especially for mas- 
sive stars where the treatment of convection, the inclu- 
sion or neglect of mass loss, the explosion mechanisms, 
etc could give rise to large discrepancies in the final yields 
of several elements. 

Most of the available models use metallicity- 
independent yields and/or concentrate on a limited num- 
ber of chemical species. TWW1995 were the first to con- 



sider the evolution of the complete ensemble of light and 
intermediate-mass elements (from H to Zn) by including 
in a simple infall model for the galactic disk the yields of 
massive stars with metallicities between Z=0 and Z=Zq 
of Woosley & Weaver (|1995|) , and adopting for the CNO 



yields of low and intermediate mass stars the Renzini & 
Voli ( |1981|) results. Recently, G P|2000| have reanalyzed the 
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evolution of the elements from C to Zn by means of an 
infall model that treats separately the halo and the disk, 
but they only take into account stellar yields from mas- 
sive stars, assuming zero net yields for intermediate mass 
stars. 

In this paper we present new calculations of the evolu- 
tion in the solar neighborhood of all the chemical elements 
up to Zn , in the framework of a two-infall model (Chi- 



appini ct al. 1997) for the formation of the Galaxy, with 
metallicity dependent stellar yields for the whole range 
of stellar masses considered. Besides, we include the 
contributions for novae nucleosynthesis, since important 
amounts of several isotopes (mainly 7 Li, 13 C, 15 N, and 
17 0) could be produced by nova outbursts. The plan of 
the paper is as follows: In Sec. 2 we present and discuss the 
main ingredients of our model: infall assumptions, SFR, 
IMF and nucleosynthetic prescriptions. Inspired by the 
recent observations of Wakker et al. ( 1999 ), who have re- 



ported the detection of a massive cloud with a metallicity 
^0.09 times solar falling into the galactic disk, we consider 
as our standard model one in which first the halo and the 
thick disk of the Milky Way form by accretion of primor- 
dial gas, and then, in a second infall episode lasting up to 
nowadays, the thin disk assembles from slightly enriched 
extragalactic material; however, we have also considered 
a more standard model with infall of primordial gas along 
the whole galactic evolution. In Sec. 3 we confront our re- 
sults with the main observational constraints in the local 
disk (G-dwarf distribution, age-metallicity relation, solar 
abundances, supernova rates, etc.). The core of the pa- 
per is presented in Sec. 4, where we compare in detail the 
evolution of the different elements studied with the avail- 
able observational data, including the most recent obser- 
vations. Finally, in Sec. 5 we give the main conclusions of 
our work. 



2. The chemical evolution model 

We have developed a standard open chemical evolution 
model where the Milky Way builds gradually up by infall 
of primordial and, lately, slightly enriched gas. No out- 
flows are considered. The disk is divided into concentric 
independent rings 1 kpc wide, and we neglect any radial 
flow between them. Inside each zone we assume instan- 
taneous mixing so that the stellar ejecta are completely 
mixed with the interstellar medium as soon as the stars 
die; in this way, each ring is composed by an homogeneous 
mixture of gas, stars and stellar remnants, and the local in- 
terstellar medium is characterized at any time by a unique 
composition; therefore, the quantities describing the state 
of each zone, i.e. surface density of total, gas and stellar 
masses, chemical abundances, etc., are only functions of 
galactocentric radius and time. We relaxed the instanta- 
neous recycling approximation by treating in detail the 
delay in chemical enrichment due to the finite stellar life- 



times, which we adopt from the work of the Geneva group 
(Schaller et al. [199$ Charbonncl ct al. |1996Q . 

The model solves numerically the classical set of 
nonlincar-integro differential equations of galactic chem- 
ical evolution (Tinsley 198C| ; Pagel 1997 ). In the solar 
ring we follow the evolution of 76 isotopes from hydro- 
gen to zinc. To keep handy the computation we use, as 
Gaussian quadrature summation for the 



in TWW1995 



mass integrals (Press et al. 1992), and a Cash-Karp step- 
per method (Press & Teukolsky 1992) for the explicit time 
integrations. In the following subsections we discuss the 
main ingredients of the model. 

2.1. Infall 

Closed box models of galactic evolution face the so called 
"G-dwarf problem" , the formation of too many long lived 
stars at low metallicities. The most common way out of 
this problem is to turn to open models which consider 
that the Galaxy forms by continuous infall of extragalac- 
tic material. In fact, there are observational indications 
of current infall onto the Galactic disk from external re- 
gions in the form of High Velocity Clouds moving towards 
the Galactic disk, as first suggested by Larson ( 1972 ). Al- 
though the interpretation of such clouds as gas of extra- 
galactic origin has been a matter of debate, the observa- 



tion of both High (Mirabel |1981|) and Very High Velocity 
Clouds (Mirabel & Morras 1984) seemed to confirm the 
idea of current infall with rates of the order of 1 Mq yr . 
Recently Wakker et al. ( 1999| ) have reported the detection 
of a massive (~10 7 Mq ) cloud falling into the disk of the 
Milky Way, with a metallicity 0.09 times solar, and the 
authors give strong arguments for the extragalactic origin 
of the cloud. 

Different types of infall have been explored in several 
models, but not all of them can solve the G-dwarf problem. 
For instance, constant infall rates just balancing star for- 
mation produce too many stars at high metallicity. Most 
of the "successful" models use simple exponentially de- 



creasing infall rates (T WW 1995; Prantzos & Aubert 1995 



Thomas et al. 1998) for the assembling of the galactic disk, 



with timescales of the order of 3-4 Gyr at the solar ring, 
but recent determinations of the meta llicity distribution 
of disk sta rs by W yse & Gilmore ( 1995 ) and Rocha-Pinto 
& Maciel ( 1996 ) require, to be reproduced by the mod- 
els, longer timescales, with typical values o f 7 Gy r at solar 
galacto centri c distances (Chiappini etal. 1997 ; Prantzos 
& Silk [1998[ Boissier & Prantzos |l999| ). It is worth to 
mention that some chemodynamical models also find long 
times cales for the formation of the disk (Samland et al. 
1997b . 



Recently Chiappini et al. (1997), and lately Chang et 
al. ( |1999| ), have included implicitly the halo phase evo- 
lution trough models that assume two subsequent infall 
episodes. In the first episode the halo and the thick galac- 
tic disk form in a very short time (tt < 1 Gyr). At the 
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end of this phase, the thin disk begins to form in a second 
infall episode characterized by much longer timescales. In 
this way, the material that eventually forms the disk has 
an extragalactic origin. 



In the present work, following Chiappini ct al. (1997), 
we adopt the two-infall exponentially decreasing model. 
Therefore, the time evolution of the total surface mass 
density in the solar ring is given by 

where tt and td are respectively the timescales for the 
halo-thick disk and thin disk phases, and t max is the time 
of maximum mass accretion onto the thin disk, which cor- 
responds to the end of the halo-thick disk phase. We set 
rr=l Gyr, and the same value for t max , and we take a 
value of 7 Gyr for T£>, the timescale for the thin disk 
formation at the position of the Sun (we adopt a solar 
galactocentric distance iq = 8.5 kpc). 

The coefficients A(rQ) and B(rQ) are fixed by impos- 
ing that the current total surface mass densities of the 
thick and thin disks in the solar neighborhood are well 
reproduced by the model 

VT(rQ,t G ) 



B(r Q ) 



tt(1 - e-te/i-T) 
[a(r Q ,t G ) -a T (rQ,t G )} 

T D [\ - e -(*G-*ma*)/r D ] 



where a(rQ,t G ) and o~T{rQ,t G ) are, respectively, the lo- 
cal surface mass densities of total and thick disk at the 
present time, t G , that we take as 13 Gyr. For a(r Q,t G ) 
we adopt a value of 54 Mq pc -2 (Rana [1991| ; Sacket |1997j ) . 
The local surface mass density of the thick disk is, in fact, 
a parameter in this kind of models. Observational esti- 
mates can be obtained from studies of the density ratio 
of the thick and thin disks, together with val ues o f their 
respective sc ale h eights (Kuijke n fc Gilmore 198S; Rci d 
& Majewski |1993| ; Robin et al. |l996| ; Buser et al. |l998| ). 
Unfortunately, the results obtained in those studies span 
a rather wide range, from 4.85 to 14.1 Mq pc~ 2 . In our 
calculations we adopt a value of 10 Mq pc~ 2 



Based on the observations of Wakker et al. (1999), we 
begin accreting primordial material during the halo-thick 
phase, and then, when the thin disk initiates its forma- 
tion, we assume that the infalling material has already 
being slightly enriched, with a typical metallicity of 0.1 
Zq in solar proportions. By comparison, we have also cal- 
culated models that only accrete primordial matter during 
the whole evolution. Even if there are not substantial dif- 
ferences between the results we obtain in both types of 
models, in agreement with Tosi ( 198§D , who showed that 
as long as the metallicity of the accreted material remains 
below 0.1 Zq there are little changes in the model results, 
we obtain slightly better agreement with the observational 
constraints in the solar region for the model that incorpo- 
rates enriched infall during the assembling of the thin disk. 



2.2. Star formation rate 

In view of the difficulties to understand the rather com- 
plicated process of star formation, standard models for 
the chemical evolution of the Galaxy adopt different ana- 
lytical prescriptions for the star formation rate (SFR) in 
terms of intrinsic parameters of spiral galaxies. The sim- 
plest and still commonly used law for star formation is the 
Schmidt ( |1959| ) law: <I> oc cr^, proportional to some power 
(between 1 and 2) of the surface gas density. Observations 
by Kennicutt ( |1998| ) of the correlation between average 
SFR and surface gas densities (total, atomic plus molecu- 
lar) in spiral and starburst galaxies point to a value of the 
exponent in the Schmidt law of fc^l.5. Besides this global 
behaviour, the SFR has to show also dependency on the 
local environment, which in turn is typically a function of 
the galactocentric distance. 

The precise radial dependence of the SFR varies ac- 
cording to the scenario considered for star formation. 
In theories that describe star formation as a local, self- 
regulating process through the balance between the grav- 
itational settling of the gas onto the disk, that enhances 
star formation, and the energy injected back into the in- 
terstellar medium by massive young stars under the form 
of winds and supernova explosions, which heat and expand 
the gas , redu cing the process of star formation (Talbot & 
Arnett 1975 ), the SFR is a function of the local gravita- 
tional potential and, therefore, of the tot al sur face mass 
density, a. In the original Talbot & Arnett (1975) formula- 



tion, \& oc a "a- We notice that some chemodynamical 



models (Burkert et al. 1992| ) find a similar dependence of 
the SFR on the total surface mass density. 

Observational support for such star formation law has 



been obtained by Dopita & Ryder (1994), who showed the 
existence in spiral disks of an empirical link between the 
H a emission, tracing current star formation, and the I- 
band surface brightness, a measure of the contribution of 
the old stellar component and, therefore, of the total sur- 
face mass density. This observed relation is well fitted by 



a SFR law very similar to that of Talbot & Arnet (1975): 
^ cx cr"a™, with n = 1/3 and m = 5/3. Moreover, self- 
regulated star formation can accommodate the observed 
correlation between surface brightness and o xygen abun- 
dance in late spiral disks (Edmunds & Pagel 1984 ; Ryder 
1995b . 



We incorporate in our numerical code this law for the 
star formation as 



a n (r,t)a 7 a n (r,t) 
v \ - 9 ' M pc" 2 Gyr" 1 



a n+m - 1 (r Q ,t) 



where the denominator is introduced as a normalization 
factor in order to express the efficiency coefficient v in 
Gyr^ 1 . Here we adopt v =1.2 Gyr -1 . 
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2.3. Initial mass function 

The initial mass function (IMF) constitutes another basic 
input in models of chemical evolution, since it determines 
in which proportions stars of different masses enter into 
play; that, in turn, fixes the averaged stellar yields and 
remnants masses for each generation of stars. 

In general, the assumed form of the IMF is a declining 
function of mass in terms of a power law: $(M) oc M~ x , 
constant in space and time. Salpeter's ( 1955| ) version, with 
x = 1.35 for the whole mass range has been commonly 
used. Ho wever , recent studies revi ewed by Scalo ( 1998Q , 
Kroupa (ll998p and Meyer et al. (|2000|) indicate, within 
still rather large uncertainties, that observations are con- 
sistent with an IMF almost flat at low masses, how flat 
s till b eing a matter of debate (for instance, Reid & Gizis 
(|1997|) find a uniqu e slope of x = 0.05 for O.1<M/M <1, 
but Kroupa (|l998| ) claims that such slope is only appro- 
priate for M/Mq<0.5), and declining as a power law with 
a slope similar to the Salpeter's one above 1 Mq, although 
there are again discrepancies on the actual value of x. For 
instance, Massey et al. (1995) find slopes ie~1-1.5 in OB 



associations, while for massive field stars Massey ( 1998 ) 
obtains very steep IMFs, with values of x~3-4. 

Another aspect of the IMF that has not yet being sat- 
isfactorily settled, nor observationally nor theoretically, 
is its time behavior. There seems to be a tendency be- 
tween observational researchers to favor no variations of 



the IMF, although in a recent review Scalo ( 1998 ) argues 
against its universality in space and time. Besides, an IMF 
producing more massive stars in the early Galaxy has been 
invoked as a solution for the G-dwarf problem. Neverthe- 
less, in a recent paper, Chiappini et al. ( 2000] ) have inves- 
tigated the effects on galactic chemical evolution of several 
time dependent IMFs and conclude that the combination 
of infall and a constant IMF is still the best choice to re- 
produce the observational constraints in the Milky Way. 

In this pa per w e adopt the constant IMF version of 
Kroupa et al. ( 1993 ) which consists of a three slope power 
law. In the range of very low masses the slope is quite flat, 
x=0.3 for M<0.5 Mq, it steepens to x=1.2 in the interval 
0.5 M<1 Mq, and in the high mas s regi me M>1 Mq the 
slope agrees with the one by Scalo (1986) x=1.7. As usual, 
we normalize to unity this IMF between a minimum stellar 
mass of 0.08 Mq, the H-burning limit, and a maximum of 
100 M . 

2.4- Nucleosynthesis prescriptions 

The chemical yields synthesized by stars in different stel- 
lar mass ranges are a key ingredient when trying to un- 
derstand the chemical evolution of galaxies, in particu- 
lar those corresponding to "massive stars" , stars that end 
their lives in Type II supernova explosions, since they are 
the main responsibles for the enrichment of the Universe in 
intermediate and heavy nuclei, through both hydrostatic 



burning (elements up to calcium) and explosive nucleosyn- 
thesis (iron peak elements, to which Type la supernovae 
are also important contributors). Even if the uncertain- 
ties in the theory of stellar evolution are far less important 
than those affecting the SFR or the IMF, there still remain 
serious discrepancies in the literature about the treatment 
of crucial aspects of stellar evolution. For instance, pre- 
supernova configurations are affected by the precise for- 
mulation of convection, scmi-convcction and overshoot- 
ing, the adopted nuclear reaction rates (especially for the 
12 C (a, 7) 16 O reaction), the inclusion or not of mass loss 
and rotation, etc. A similar situation holds for the results 
of core collapse supernova explosions. In fact, in current 
calculations the explosion itself is induced rather arbitrar- 
ily and the explosion energy is fixed by imposing a given 
value for the final kinetic energy. Also, the ejected mass 
crucially depends on the degree of "fall-back" , and thus 
on the precise details of the explosion. We thus see that 
the published stellar yields are yet plagued with severe 
uncertainties. 

Although the evolution during the thin disk epoch will 
be hardly affected by the slight metallicity dependence 
shown in the published stellar yields, the early halo-thick 
disk phase could be strongly influenced by the material 
ejected by stars of different metallicities. Hence, we will 
consider metallicity dependent yields for the whole stellar 
mass range. 



2.4.1. Massive stars: 8 M <M<100 M 



*0 



We assume that stars in this mass interval end their lifes 
as core collapse supernovae. The two major sources in this 
mass range, widely used in ch emica l evolution mod els, ar e 
those of Woosley & Weaver ( |l995|) (hereafter WW |l995|) , 
that calculated full stellar models without mass loss nor 
rotation, for stars with masses comprised between 12 
Mq and 40 Mq and with different initial metallicities 
(Z/Z =O, 10 " 4 , IP " 2 , lO" 1 and 1), and those of Thiele- 
mann et al. (1996), who calculated the evolution of He 



cores corresponding to stars of solar initial metallicity 
up to 7 Mq . There are noticeable differences between 
the WW|1995| results for Z=0 and Z^0, while the yields 



and stellar remnant masses for nonzero initial metallici- 
ties give more similar results. Solar composition models 
of WW199£ and Thielemann et al. (1996) give almost the 
same values for the yields of the CNO isotopes, but non- 
negligible differences exist for other important elements. 
For instance, the magnesium yield is systematically lower 



in WW 199!: as a consequence of the different treatments 



of convection. In the case of the iron yield, WW 1995 gives 



much higher values than Thielemann et al. (1996) for stars 
with masses below ^35 Mq, since it is very sensitive to 
the mass cut that separates the ejecta and the material 
that falls back onto the compact remnant which, as men- 
tioned above, depends in turn on the detailed explosion 
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mechanism (see Thomas et al. 1998, and Chiappini et al. Mq, besides of taking into account convective overshoot- 



1999 for a complete discussion) 



Lately, Limongi et al. ( 2000 ), from full stellar models, 
again without mass loss nor rotation, have also presented 
yields of massive stars in the range 13-25 Mq for three 
metallicities (Z=0, 10~ 3 , and 0.02). As sh own b y GP|200C, 



the main difference with respect to WW 1995 is a more 
marked odd-even effect in the case of Limongi et al. ( [2000 ) 
yields, which translates into systematically lower yields of 
odd Z elements. 

In t his pa per we adopt the metallicity-dependent yields 
of WW 1995] for massive stars between 8 and 100 Mq. The 
reasons for this choice are twofold: they are the results of 
full evolutionary calculations and they consider a wider 
interval of masses and metallicities as compared with the 
work of Limongi et al. ( |2000|) . In particular, we take their 
models A up to 25 Mq, and their models B for 30, 35 and 
40 Mq (the explosion energy in models B is higher than 
in models A). Those yields have been extrapolated up to 
100 Mq, even though this has a minor effect on the final 
results since the Kroupa et al. ( 1993] ) IMF produces very 
few stars more massive than 40 Mq. 



TWW |1995| suggested that a better agreement with 
most of the observed evolution of the abundances is ob- 



tained if the WW |1995| iron yields are reduced by a factor 
of two. Observational estimates of the iron synthesized in 
SN 1987A (the explosion of a 20 Mq belonging to a stel- 
lar system, LMC, with an estimated Z~0.1 Zq) and of SN 
1993J (whose progenitor was a 14 Mq star in the galaxy 
M81, where Z ~Zq) also points to such an overestimate 
(Thomas et al 1998). Therefore, in view of the uncertain- 
ties arising from the flaws in the explosion calculations, 



we reduce by a factor 2 the nominal yields of WW 1995 



2.4.2. Low and Intermediate Mass Stars: M<8 Mq 

As it is well known, single stars in this mass range pollute 
the interstellar medium through moderate stellar winds 
and planetary nebula ejection, ending their lifes as white 
dwarfs. These stars make important contributions to the 
He, C and N galactic contents. Here again important 
uncertainties on the final yields exist that result from 
the treatment of mass loss, convection, evolution on the 
asymptotic giant branch, etc. 



The yields of Renzini & Voli (1981) have been exten- 
sively used in chemical evolution models. However, in the 
last years several groups have published new evolutionary 
calculations for stars belonging to this mass interval and 
for different initial metallicities (Marigo et al. 1996; van 



ing, which is not considered in WW 1995 



2.4.3. Type la Supernovae 

There is a general agreement that Type la supernovae are 
important contributors to iron and other iron peak ele- 
ments (in particular 58 Ni and 54 Cr) in the late disk evo- 
lution, up to the point that they are r esponsible for ^2/3 
of the total iron contents (T WW 1995). We adopt the con- 
ventional view that Type la supernovae are the carbon de- 
flagration of massive C-0 white dwarfs in binary systems. 
The contribution of Type la supernovae to the chemical 
enrichment of the galaxy, specially the iron content, has 
been calcul ated f ollowing the prescriptions of Matteucci 
& Greggio ( 1986 ), where Type la supernovae come from 
binary systems with a minimum mass of 3 Mq , in order to 
ensure that the accreting white dwarf eventually reaches 
the Chandrasekhar mass, and a maximum mass of 16 Mq, 
if we assume that C-0 white dwarfs come from primaries 
up to 8 Mq. Thus, the amount of isotope i due to this 
source is given by 



C J $(M B ) 



0.5 



f( f i)y(r,t-T fl .M B )jj-dn 



dM B 



The parameter C, which determines the fraction of 
suitable binary systems that actually produce Type la SN, 
is set equal to 0.04 in order to fulfill the restrictions im- 
posed by observational data on supernovae frequency. The 
chemical composition of the ejecta are taken from Thiele- 



den Hoek & Groenewegen 1997). In this work the contri- 
bution to the galactic enrichment by these stars has been 

taken from van den Hoek & Groenewegen (1997), who have masses between M™ =0.5 Mq and M^ 1 =1.5 M 



mann et al. (1993) for the well known model W7 (Nomoto 
et al. 1984 ), and we also include the contribution to the 
enrichment of the interstellar gas by the secondary star. 



2.4.4. Classical Novae 

Although classical novae have little or no influence on the 
evolution of the abundances of most of the chemical el- 
ements included in our calculations, according to Jos & 
Hernanz ( 1998) they can be significant contributors to cer- 
tain isotopes that are strongly overproduced in nova ex- 
plosions relative to the solar system abundances (in partic- 
ular those with overproduction factors larger than 1000). 
7 Li, 13 C, 15 N and 17 are some of the isotopes whose 
abundances could be affected by novae. 

In order to obtain the nova rates we assume that novae 
take place in binary systems made of a white dwarf com- 
ing from stars with main sequence masses in the interval 

My 



1 Mq to Mf =8 Mq, while the secondary stars 



give yields for stars with masses comprised between 0.9 
and 8 Mq and consid er five initial metallicities, since the 
work of Marigo et al. ( 1996 ) is limited to stars with M<4 



©■ 

The rate of those explosions can be estimated by a proce- 
dure similar to that used by Matteucci & Greggio ( 1986 ) 
to calculate the rate of SN la: 
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dRou 



tburts 



(It 



= D 



<f>(M B ) 
M B 



dMr 



M(t)+M r 2 

f(jjL)ip{t - r Mfl(1 _ Al ) - t coo i)dn 



where Mb is the total mass of the system and 
M B > 2M 2 M 



Mm 




M B < 2M$* 



( m b -m™ Mi\ m m m 

\ M B ' M B J m B > m l 



M, 

m: 
m' b 
0.5 



MJT- 



> M B > 2M 2 m 

> M B 



In this expression, the value of the free parameter D, 
which plays a role similar to the parameter C in the rate of 
Type la supernovae, is obtained by imposing that the for- 
mula reproduces the observed outburst rate at the present 
age of the Galaxy, of the order of ^40 yr~ 1 (Hatano et al. 
1997), and t coo ;, the cooling time for a white dwarf before 



it can produce the first nova outburst, is taken as 1 Gyr . 

The yields for the classical nova outburst are an av- 
erage of those given by Jos & Hernanz ( 199§| ). We have 
considered that 30% of the nova outbursts come from ONe 
white dwarf progenitors, while 70% come from CO ones. 



3. Model results versus observational constraints 

As we have discussed in the previous section, there are 
large uncertainties in the ingredients of standard chemical 
evolution models which translate into a series of more or 
less free parameters, e.g. infall timescales, efficiency and 
exponent of the SFR, slope of the IMF, etc. Fortunately, 
for the solar neighborhood, the number of observational 
constraints is large enough to restrict the space of varia- 



tion of these parameters (Boissier & Prantzos 1999). Re- 
sults for the whole galactic disk will be the subject of 
another paper. 

3.1. Surface densities. SFR and supernova rates 

First of all, any successful model has to reproduce the 
present-day surface density of total mass, gas, star and 
stellar remnants, the current rate of star formation, infall 
and supernovae (both Type la and Type II). In Table 
we show our results for the present values (t<3=13 Gyr) of 
those quantities compared with the current observational 
values for the solar vicinity. The time evolution of the 
surface densities of total mass (a), stars (<r*)> S as {°~g)i an d 
stellar remnants (<r rem ) are displayed in the upper panel of 
Fig. [l], while in the lower panel appear the star formation 
rate, and the Type II and Type la supernova rates. We 
see from Table [j] and Fig. [j] that the model presented in 
the previous section nicely fits the observations. 



Table 1. Calculated and observed present values for sev- 
eral relevant quantities in the solar ring. 





model 


Obser. 


Ref. 


a g (Mq pc a ) 


8.27 


7-13 


1 


a, (Mq pc" 2 ) 


39.93 


30-40 


2 


O r em (Mq pC 2 ) 


5.77 


2-4 


3 


Og/lJtot 


0.15 


0.05-0.2 




(Tin/ (M pc -2 Gyr" 1 ) 


1.38 


0.3-1.5 


4 


* (M pc" 2 Gyr" 1 ) 


2.84 


2-5 


5 


SN la per century 


0.41 


0.17-0.7 


6 


SN II per century 


1.76 


0.55-2 


6 


SN II/SN la 


4.3 


3-6 


6 


Novae (yr~ ) 


40.0 


20-60 


7 
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Fig. 1. Upper panel: time evolution of the surface densities 
of total mass (a), visible stars (c*), gas (cr 5 ), and stellar 
remnants (a rem ) for our two-infall model. Lower panel: 
evolution of the stellar formation rate, as well as the Type 
II and Type la supernova rates. 



3.2. Age-metallicity relation 

Our model also fulfills the requisite of producing an age- 
metallicity relation in agreement with the observed one. 
The age-metallicity relation (AMR) shows the evolution 
of the ratio [Fe/H], taken as a measure of the metallicity 
of the Galaxy, as a function of time. The pioneering work 
by Twarog p80| ) showed that [Fe/H] rises steeply during 
the first 2 Gyr up to a value of ~— 0.5 dex, and afterwards 
increases smoothly with time towards the solar value, al- 
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though the data showed a large dispersion, and age-bins 
and average metallicity per bin were used. Later reexami- 



nation of Twarog's data by Meussinger et al (1991) found 
a similar AMR. More recently, Edvardsson et al. ( 1993Q , 
based on a different sample of nearby F and G stars, ob- 
tained for binned data a good agreement with previous 
results, but they claim that the large dispersion found in 
their unbinned data was essentially real and not obser- 
vational (but see Garne tt & Kobulnicky 2000). Finally, 



Rocha-Pinto et al. ( 2000 ) obtained a chromospheric AMR 
using a sample of 552 late-type dwarfs that basically con- 
firms the observed trends, b ut wit h almost half the scatter 
found by Edvardsson et al. (1992 ) . 

Models of chemical evolution that adopt the instan- 
taneous mixing approximation, like ours, cannot produce 
any dispersion at all, and at most they can only fit the 
mean relation. Therefore, the AMR is a less tight con- 
straint, even though the average relation must be fitted 
by successful models. 

We present in the lower panel of Fig. || our results for 
the evolution with time of [Fe/H] compared with the AMR 
by Meussinger et al. ( 199l| ), Edvarsson et al. (1993) and 
with the more recent data of Rocha-Pinto et al. ( 2000) ). 
Just for completeness, in the upper panel we also show 
the evolution of the total metallicity in terms of the solar 
value. Again our model fits well the averaged trend of 
the data, although, as stated above, this is not a really 
tight constraint, because almost any model can produce a 
reasonable AMR. 

Nearly 2 /3 of the present total iron contents are made 
by Type la supernovae, the other third coming from Type 
II supernovae. We recall here that the iron yields for Type 



II supernovae used in this work are half those of WW 1995 



in order to obtain a better fit to the evolution of the el- 
ement abundances. That is also indicated by iron abun- 
dance measurements in supernova ejecta. In any case, as 



pointed out by TWW1995 whatever contribution to the 



present galactic iron contents by Type II supernovae com- 
prised between one and two thirds is allowed by our cur- 
rent understanding of stellar physics. 

3.3. G- dwarf distribution 

The stars of spectral type G have main-sequence lifetimes 
comparable or even larger than the estimated age of the 
Galaxy. Hence, the distribution in metallicity of a com- 
plete sample of these stars in the solar neighborhood car- 
ries memory of the star formation history. This makes the 
observed G-dwarf distribution one of the more stringent 
constraints for models of galactic chemical evolution. 

As it is well known, the paucity of G-dwarf stars at 
low metallicities, referred to as the "G-dwarf problem", 
cannot be explained by simple, closed box models. Sev- 
eral solutions have been proposed to solve this problem, 
but open models with progressive infall of primordial or 
slightly enriched material with long timescales for the disk 
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Fig. 2. Upper panel: evolution with time of the total 
metallicity Z normalized to the solar value Zq. Lower 
panel: time evolution of the ratio [Fe/H] for contributions 
from Type II supernovae (thin line) and for Type II plus 
Type la supernovae (thick line). 



formation, like ours, are still the best option (Chiappini 
et al. 



1997), besides of being compatible with dynamical 



simulations of the formation of galactic disks (Burkert et 
al. 199S ) , and with the observation of infall of High and 



Very High Velocity Clouds. 

Fig. H compares the predicted G-dwarf distribution 
produced by our model with the most recent observational 
data from Wyse fc Gil more ( 1995j ) (left panel), Rocha- 
Pinto & Maciel ( 1996|) (central panel) and J0rgensen 
( 2000| ) (right panel). The direct results are displayed by 
thin dashed lines, while the thick solid ones show the con- 
volution with a Gaussian with a dispersion of 0.15 in order 
to simulate observational and intrinsic scatter. Very good 
agreement is reached, specially for the rising part of the 
distribution and the location and height of the peak when 
comparing with the two first sets of data. In the case of 
the distribution from J0rgensen ( |2000| ), the fit is still rea- 
sonable, the peak is again well reproduced, but the model 
shows a low metallicity tail (around [Fe/H]^ —0.5) that is 
higher than the data. The comparison between theoretical 
results and observations could be improved by relaxing the 
instantaneous mixing approximation, introducing thus in 
the models the intrinsic scatter. Our results confirm once 
more that infall models with long timescales for the as- 
sembling of the disk are capable of solving the G-dwarf 
problem. 
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Fig. 4. G-dwarf distributions obtained for different types 
of infall. The obs ervati onal data correspond to Rocha- 
Pinto and Maciel fll996|). 



As indicated above, our model uses a double exponen- 
tial infall in order to treat as separate entities the halo- 
thick disk and the thin disk, and we consider accretion of 



enriched material (Z=0.1Zq) during the thin disk phase. 
We have also calculated models that incorporate primor- 
dial material both in the halo-thick disk and in the thin 
disk phases, as well as models that only consider simple ex- 
ponential infall of primordial and primordial plus enriched 
material which, if not appropriate for the halo and thick 
disk, are relevant when making comparisons for the thin 
disk properties. In all those models we have adopted the 
same timescale for the galactic disk, i.e., 7 Gyr. We obtain 
in all cases very similar results for the final characteristics 
of the solar ring, but slight differences appear for the con- 
volved G-dwarf distributions, which are displayed in Fig. 
^, where for clarity only the observations by Rocha-Pinto 
& Maciel ( 199(f ) are shown. One and two infall models 
with the same composition of the incorporated matter are 
almost indistingishable, but primordial composition mod- 
els show a moderate excess of stars in the low metallicity 
tail and lower maxima than those obtained in enriched 
models. 



3.4- Solar abundances 

We present in Fig. [| the calculated mass fractions of the 
elements (upper panel) and all their stable isotopes (lower 
panel) studied in this work (except LiBeB isotopes, whose 
abundances have important contributions from galactic 
cosmic ray nucleosynthesis, not included in the present 
calculations) in the interstellar gas at the time of the so- 
lar birth, 4.5 Gyr ago, compared with the observed solar 
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Fig. 5. Mass fractions of the calculated elements (upper panel), from hydrogen to zinc (lithium, beryllium and boron 
are not shown), and of all their stables isotop es (lo wer panel) in the interstellar gas at the time of solar birth, relative to 
the solar abundances of Anders & Grevesse ( 1989 ). The isotopes of the same element are connected by solid lines. The 
dashed horizontal line would correspond to a perfect agreement, while the two solid lines correspond to a discrepancy 
by a factor of 2. We take as successful those elements and isotopes whose calculated abundance fall within this factor 
of 2 of their solar value. The three open circles in the lower panel represent the relative abundances of 13 C, 15 N and 
17 when nova yields are excluded from the calculations. 



values (Anders & Grevesse |1989[ ), considered as represen- 
tative of the composition of the interstellar gas at that 
epoch. We must recall here that this last assumption is far 
from being definitely settled, since measured abundances 
in regions of recent star formation, like the Orion neb- 
ula, show metallicities lower than solar (see, for instance, 
Cardelli & Federmann |l997[ ). 

Due to the uncertainties existing in the observational 
values, we consider as good agreement abundances within 
a factor of 2 of their solar value. As it can be seen in Fig. 
||, all the elemental abundances fulfill this criterion, nickel 
excepted, which is clearly oversolar, and the model makes 



the correct solar metallicity. The isotopic abundances are 
also nicely adjusted. The spread for isotopes below calcium 
is lower than above it, which reflects the uncertainties in 
the present modeling of Type II supernova explosions. 

The production factors for hydrogen, helium and the 
most important CNO isotopes are nearly equal to unity. 
In particular, the solar 3 He abundance is almost exactly 
reproduced by our models without taking into account any 
partial destruction by extra-mixing on the RGB (see, how- 
ever, Romano et al. 2000). The inclusion of novae wors- 
ens our results for 13 C and 17 0, which are then a little 
bit overabundant. On its turn, 15 N is a factor of 3 be- 
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low its solar value when only nucleosynthesis in massive 
stars is considered, and contribution by novae is manda- 
tory to render its abundance compatible with the solar 
abundance. 

A few isotopes fall out of the acceptable range. 48 Ca, 
47 Ti, 50 Ti and 64 Zn, are underproduced (especially 48 Ca). 
In the case of the neutron-rich isotopes 48 Ca and 50 Ti 
there could be a substantial contribution from Type la 
supernovae resulting from initially slow deflagrations at 
high density (Woosley & Eastman 1994). 47 Ti (and, in 
lesser extent, 44 Ca and 51 V) have been traditionally prob- 
lematic; these isotopes could have noticeable contributions 
from sub-Chandrasekhar models of SNIa (Shigeyama et al. 



1992; Woosley & Weaver 1994). Additional contributions 



to the abund ance of 64 Zn could be obtained, as pointed 
out by T WW 1995, from weak s-process in massive stars 
or classical s-process in low-mass AGB stars. 

On the other hand, there is a small overproduction of 
nickel, due to the large yields of 58 Ni 58 Ni in the adopted 
model for SNIa and the high production of 62 Ni 62 Ni in 
current models of massive stars explosions which reflects 
the inaccuracies involved in present supernova modeling. 



As noted in previous works (TWW1995; GP2000), the 



fact that the solar abundances of all isotopes up to the 
iron peak, that cover a range of almost 9 orders of mag- 
nitude, are so precisely reproduced by current models of 
chemical evolution indicates that the present understand- 
ing of stellar nucleosynthesis is basically correct, at least 
to first order. 

4. Abundance ratios evolution 

We have calculated the evolution in time of all the stable 
isotopes of thirty chemical elements from H to Zn with 
the nucleosynthetic prescriptions detailed above. Our re- 
sults are displayed in Figs. ^ to [l5[ where we compare the 
predicted ratios [X/Fe] as functions of [Fe/H] with obser- 
vational data of spectroscopic abundances in nearby stars. 
Data sources are detailed in the legends of the figures. As 
the data normally refer to elemental abundances, for each 
element we sum over the calculated isotopic abundances. 

Before going on with the presentation of our results 
we wish to recall again that following the suggestion of 



TWW1995 we have reduce the iron yields of WW 1995 



by a factor of 2. The reduced iron yields, well within the 
current inaccuracies in supernova models, not only give a 
better agreement with the observed elemental abundances 
but are also consistent with the amount of iron ejected in 
recently observed Type II supernovae . 

We consider that the calculated abundance evolution 
is only significant for [Fe/H]>— 3, a value that our model 
reaches at t~30 Myr, which approximately corresponds to 
the lifetime of a 8-9 Mq star. This means that only after 
reaching such global metallicity the yields of massive stars 
are fully averaged over the IMF, therefore reducing the 
possible uncertainties in the yields of individual stars. 




-2.0 -1.0 
[Fe/H] 

Fig. 6. Evolution of abundances ratios [X/Fe] as a func- 
tion of [Fc/H] for the CNO elements. Upper panel: Carbon; 
central panel: Nitrogen; lower panel: Oxygen 



4-1- The CNO elements 

We present in Fig. || the results obtained in our calcula- 
tions for the evolution of CNO elemental abundances with 
respect to iron. 

Carbon and iron are primary elements and are pro- 
duced by very different processes: carbon is made by the 
triple-a process in hydrostatic helium burning in massive 
and, mostly, in intermediate and low-mass stars, while Fe 
is produced by explosive burning in Type II and Type 
la supernovae. Therefore, we expect that the ratio [C/Fe] 
remains nearly constant as a function of [Fe/H]. Observa- 
tions of carbon abundances in halo and disk dwarfs (field 
giants are not reliable indicators of ab initio carbon abun- 
dance due to the effects of the first dredge-up) do actually 
show that [C/Fe] is almost constant in time, with a value 
approximately solar (Laird 1985; Tomkin et al. |1986; 
bon et al. [1987 ; 



Car- 



1995| ), although there 



McWilliam et al. 
is an important dispersion in the data at all metallicities. 
Those surveys also seem to indicate some trend towards 
slightly higher values of [C/Fe] at low [Fe/H], as confirmed 
to first order by Wheeler et al. (1989) after reanalyzing the 
data available at that epoch, although observational un- 
certainties counsel to take with some caution this upturn 
at low metallicity. 

Our results show that the synthesis of carbon and iron 
by massive stars suffices to explain the observations in 
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halo stars. The [C/Fe] ratio shows a slow decline towards 
the solar value during the halo-thick disk epoch, which in- 
dicates that intermediate-mass stars, whose lifetimes per- 
mit them to evolve during the late halo phase, do not give 
very important contributions to the carbon contents. The 
bump around [Fe/H]~— 1, at the beginning of the thin disk 
phase, results from the contribution of intermediate and, 
specially, low-mass stars (M < 2 Mq), which eject impor- 
tant amounts of carbon but no iron. But the iron ejected 
by Type la supernovae, which start to contribute around 
the same epoch, compensates and finally overwhelms the 
role of the low mass stars. 

The synthesis of 13 C in intermediate and low-mass 
stars is capable of building up the solar contents. The 
inclusion of novae produces a small overabundance of 13 C 
at the time of the solar birth. 

The evolution of the ratio [N/Fe] as a function of 
[Fe/H] is displayed in the second graph of Fig. |[ The three 
observational surveys shown in the figure present a large 
scatter and there are no observations below [Fe/H] ^—2.5. 
The data seem to indicate that [N/Fe] remains more or 
less constant over the range of metallicities studied, but 
besides the large dispersion, there are also uncertainties 
on the value of the constant, although the data are com- 
patible with a solar value. Our model predicts a rapid in- 
crease of [N/Fe] at low metallicity, which show that metal 
poor massive stars do not produce primary nitrogen. The 
ratio [N/Fe] steadily increases up to [Fe/H]~ —1 due to 
the progressive contribution of mostly secondary nitrogen 
ejected by intermediate and low-mass stars, until the iron 
production by Type la supernovae compensates the ejecta 
of intermediate stars, flattening the evolution of [N/Fe]. 

T WW 1995 obtained primary nitrogen production by 
metal poor massive stars with masses above 30 Mq by 
considering enlarged convective overshooting, but still 
within the theoretical uncertainties, due to the violent ef- 
fects of the penetration of the helium convective burning 
shell into the hydrogen shell, although this effect does not 
occur in more metal-rich massive stars. This could be a 
promising mechanism to explain the abundance pattern of 
damped Lyman-a systems (Matteucci et al. 1997). Syn- 



thesis of primary nitrogen due to the injection of protons 
into helium burning zones could also res ult in stellar mod- 
els that inclu de rotation (Heger et al. 1999 ; Maeder & 
Meynct 2000 ). Massive stars produce 15 N but the domi- 
nant contribution in our model comes from novae. 

Oxygen is exclusively produced by massive stars and 
dominates their ejecta. Until recent years the observa- 
tional status of the oxygen abundance in dwarfs and G and 
K giants showed a nearly constant value of [O/Fe]~0.5 for 
[Fe/H]< — 1, and a gradual decline in the disk. In fact, this 
is the canonical behavior of the so-called a-elements (O, 
Mg, Si, S Ca, Ti): almost flat evolution in the halo, and a 
gradual decline in the disk due to the iron production by 
Typ e la s upernovae. However, re cent d ata by Israelian et 
al. ( 1988 ) and Boesgalard et al. ( 199E ), which agree with 
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Fig. 7. Upper panel: evolution of the ratio [F/Fe] as a 
function of [Fe/H]. Lower panel: evolution of the ratio 
[F/O] as a function of [Fe/H], compared with observa- 
tional points from Jorissen et al. (1992). 



earlier results from Abaia & Rebolo (1989), are in contra- 
diction with the idea of an oxygen plateau, since they find 
that [O/Fe] increases between [Fc/H]=— 1 and -3, from 0.6 
to 1, which means a slope in this range of -0.31±0.11. 

Our calculations show that [O /Fe] begins with a value 
of the order of 1 at very low metallicity, and it slowly de- 
clines during the thick disk phase, in part because massive 
stars of different masses and initial metal contents produce 
different O/Fe ratios, and also because Type la super- 
novae begin to inject iron already in the halo-thick disk 
phase, though their influence is stronger during the thin 
disk epoch, when Type la supernovae reach their highest 
rate, as it is clearly shown by the steeper decline of [O /Fe] 
for [Fe/H]>-1. This behaviour is in agreement with previ- 
ous results by TWW |l995| and Chiappini et al. ( |1997|) . We 
obtain a slope of -0.28, which clearly indicates a devia- 
tion from the traditional oxygen plateau, although a little 
bit less pronounced than in the new data. We will see the 
same kind of evolution, although less marked, for the rest 
of the a-clcments. 

4-2. Fluorine 

In our calculations fluorine is synthesized (together with 
7 Li and 11 B which will not be treated here) by the so- 
called neutrino-induced nucleosynthesis in Type II super- 
nova explosions. The amount of fluorine produced by neu- 
trino spallation of 20 Ne is extremely sensitive to the neu- 
trino fluxes and spectra used in the explosion calculations, 
which render its yield very uncertain. 

The upper panel of Fig. ^ displays the evolution versus 
iron abundance of the ratio [F/Fe]. Our model predicts a 
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practically constant and almost solar [F/Fe]. In the lower 
panel of Fig. [?] we present the ratio [F/O] compared to 
the only available set of observations in the literature. 



In agreement with TWW1995, our model predicts strong 



subsolar [F/O] ratios at low metallicity, which is a result 
of the early efficient oxygen enrichment by massive stars 
and not of a deficient production of fluorine, while it fits 
the observed abundances of normal K giants. 

Most of the plotted data correspond to spectroscopi- 
cally peculiar giant stars whose abundances are affected 
by stellar evolution, although the high ratios observed in 
those stars also point to contributions other than from 



massive stars. Forestini & Charbonncl (1997) suggest fluo- 
rine production in the helium burning shell of AGB stars, 



and Meynet & Arnould (2000) contributions from WR 
stars. A reliable history of fluorine will only be settled 
when we gain new insight on the contributions from dif- 
ferent kinds of stars. 



4-3. The a-elements Magnesium, Silicon, Sulphur and 
Calcium 

Our results for the evolution of the ratios [a/Fe] for mag- 
nesium, silicon, sulphur and calcium, plotted in Fig. §L 
closely follow the general trends of the data. Chiappini 
et al. ( 1999| ) pointed out that the so-called "plateau" for 
the a-elements at low metallicities is not perfectly flat, 
but instead presents a slight slope. In fact, that is exactly 
what we find: the abundance ratios slowly increase with 
decreasing [Fe/H], in good agreement with the work of 
McWilliam et al. ( ]1995[ ) and that of Ryan et al. ( |1996| ) for 
magnesium, silicon and calcium, although the slopes are 
shallower than for oxygen. 

The results for silicon and calcium are in excellent 
agreement with the data. The evolution of magnesium and 
sulphur, however, is less satisfactory. 

In the case of magnesium, though we begin with typ- 
ical halo values, the rapid decline of [Mg/Fe] gives values 
a little lower than observed in the thick disk phase, but 
still marginally compatible with the observations. Later 
on, when the thin disk starts, the iron from Type la su- 
pernovae worses the problem and we clearly underestimate 
the [Mg/Fe] ratio in the disk. In principle the metallicity- 
independent yields for massive stars from Thielemann et 
al. ( |1996| ) or the WW 1995] yields for solar metallicity could 



solve this discrepancy (Thomas et al. |1998| Chiappini et 
al. |1999[ ), but this type of yields are not appropriate when 
studying the galactic halo. Since the magnesium yields 
of WW1995 increase with the stellar mass, the use of the 
rather steep IMF of Kroupa et al. ( |1993| ) worsen the situa- 
tion. As pointed out by GP 2000| the use of the metallicity 
dependent yields of Limongi et al. ( |2000p will probably 
solve the underproduction of magnesium but the ratios 
a/Mg will not reproduce the observations. The fact that 
current yields from massive stars cannot completely ex- 
plain the magnesium evolution could indicate the need of 
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Fig. 8. Evolution of the abundance ratios to iron of the a 
elements Mg, Si, S, and Ca as a function of [Fe/H]. 



a supplementary source for this element, either in low and 
intermediate mass stars or by enhancing the magnesium 
yield in Type la supernovae (TWW |l995| ). 

The paucity of data for sulphur at low [Fe/H] does 
not allow precise comparisons, but if the points by Fra- 
nois ( 1987 , 198S) around and below [Fe/H] ~- 1 adequately 
represent the main trend, a supplementary production of 
sulphur at low metallicity seems necessary. Given our nu- 
cleosynthetic sources, sulphur is only produced by massive 
stars, and the uncertainties in the WW 1995 yields of metal 
poor massive stars could be the cause of the low [S /Fe] at 
low [Fe/H]. It has been suggested, however, that the [S/Fe] 
values measured in the halo stars should be reduced by 0.2 
dex (Lambert 1989). If that were the case, our results for 
sulphur would nicely fit the corrected data. 



4.3.1. The Magnesium isotopes 

Unlike most of the chemical elements, in the case of mag- 
nesium the evolution of the isotopic ratios can be con- 
fronted with observational results. 
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Fig. 9. Evolution of the magnesium isotopic ratios 
25 Mg/ 24 Mg (upper panel) and 26 Mg/ 24 Mg (lower panel) 
as a function of metallicity. 



In the upper panel of Fig. |^ we show our results for 
the evolution of the ratio 25 Mg/ 24 Mg as a function of 
[Fe/H], while in the lower panel appears the evolution of 
the ratio 26 Mg/ 24 Mg. The three isotopes are produced 
during the hydrostatic evolution of massive stars. Since in 
the WW1995 yields for massive stars 24 Mg is basically a 
primary isotope, while the neutron rich isotopes 25 Mg and 



26 Mg increase with metallicity (both are affected by the 
neutron excess), it is to be expected, indeed, that both 
ratios decrease towards lower metallicities. 



Our results, similar to those found by TWW199E and 



by GP2000, follow the expected trend, but the model 



produces isotopic ratios clearly below the observations 
in the halo-thick disk phase (i.e., [Fc/H]<-1). If the 
WW 1995 yields do not underestimate the magnesium iso- 
topic yields, a supplementary source of 25 Mg and 26 Mg is 
needed, for instance through s-process in AGB stars (Ibcn 
& Renzini |l985| ). 

In the thin disk phase ([Fe/H] > -1), the model re- 
produce well the global trend of the observations, with 
both ratios systematically increasing with metallicity. The 
ratio 25 Mg/ 24 Mg reaches the solar value. However, the 
26 Mg/ 24 Mg ratio at [Fe/H]=0 is ~50% hig her th an solar, 



again in agreement with the results of GP200C who at- 



tribute the overabundance of 26 Mg at the solar birth to 



the effects of the IMF of Kroupa et al. (1993) that favor 
the production of 26 Mg in front of that of 24 Mg. 

4-4- Sodium, Aluminum and Potassium 

Sodium, aluminum and potassium are odd elements (be- 
sides, sodium and aluminum are monoisotopic) mainly 
produced by massive stars. They are subjected to the odd- 
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Fig. 10. Evolution of Na, Al and K abundance ratios to 
iron in terms of [Fe/H] as metallicity indicator. 



even effect, i.e., stellar yields increasing with metallicity. 
Therefore, their galactic abundance ratios to iron must in- 
crease with time, although the evolution of the ratio [ele- 
ment/Fe] could be affected by the iron production by Type 
la supernovae, especially in the disk. One can eliminate 
the effects of Type la supernovae on the abundance evo- 
lution of od d elem ents by plotting the ratio [element /Mg] , 
as in TWW 1995 , because magnesium is also nearly exclu- 
sively synthesized in massive stars, or by using as metal- 
licity indicator an element more reliable than iron, which 
is affected by uncertainties such as the evolution of Type 
la supernova rates or the mass cut and explosion energy in 
Type II supernovae, as done by GP 2000 who use calcium 
instead of iron as metallicity indicator. 

We display in Fig. |l^ the calculated [Na/Fe] (upper 
panel) and [Al/Fe] (central panel). In Fig. |ll| we also show 
the ratios [Na/Ca] and [Al/Ca] using [Ca/H] as a measure 
of the metallicity. It is worth to stress that the evolution 
of the calcium abundance obtained in our model agrees 
very well with the data (see Fig. ||). 

The observed sodium to iron ratios display a plateau 
in the disk. The situation for the halo stars is less clear. 
The rather old data by Peterson ( |198l[ ) indicate a rapid 
decrease from [Fe/H]= —1 to [Fe/H]= —2. Recent ob- 
servations by Stephens ( 1999| ) show the same trend, but 
most other observations point to a constant [Na/Fe]~0, 
although with large dispersion. We obtain an almost con- 
stant, but lower than the data, abundance up to [Fe/H]~- 
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Fig. 11. Evolution of the ratios [Na/Ca] and [Al/Ca] in 
terms of [Ca/H] as metallicity indicator. 
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2.5, and a steep increase at the beginning of the thin disk 
formation that ends in a plateau lower than observed. 

The calculated evolution of [Al/Fe] begins almost flat 
and is clearly higher than the data from halo stars. Later, 
it slowly increases at the thick-thin disk transition, but, 
as with Na, remains below the disk star observations. 

The abundance ratio to calcium of both elements 
shows now the odd-even effect (see Figjll]). We still have 
low [Na/Ca] in the complete range of [Ca/H], but the alu- 
minum ratio shows a better agreement with the observa- 
tions. We remind that our models underproduce the solar 
sodium by a factor of ~1.4, and aluminum by a factor of 
-1.2. 

As pointed out by TW W[1995 
intermediate-mass stars could 
and aluminum through the Ne-Na 
improving the fit to the data. 

The potassium abundance is dominated by the odd 
isotope 39 K, and it is again mainly produced in mas- 
sive stars through hydrostatic oxygen burning, so that 
it should show the odd-even effect. Surprisingly, the ob- 
served abundances display the opposite evolution: the few 
halo observations give values of [K/Fe]~0.5 with a wide 
scatter, decreasing to solar values in the disk. Our results, 
presented in Fig.jlO], show, however, a behavior similar to 
other odd elements. We obtain slightly lower than solar 
values in the halo-thick disk, and a gradual rise during 
the thin disk phase towards final solar values, in agree- 



and GP|200C, 
produce some sodium 
and Mg-Al cycles, 



ment with the reduced iron yield evolution of TWW1995 



4-5. Neon, Phosphorus, Chlorine and Argon 

Although there are not observational data to compare 
with, we plot in Fig. h2 the calculated evolution of the 
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Fig. 12. Evolution of neon, phosphorus, chlorine and ar- 
gon abundance ratios to iron, as a function of [Fe/H]. 



noble gases neon and argon, as well as those of phospho- 
rus and chlorine, just for completeness. 

Neon abundance is dominated by 20 Ne, and the three 
argon isotopes are also of even Z. Therefore, both elements 
should then behave as a-elements. The calculations con- 
firm this theoretical expectation. The evolution of neon 
in the halo-thick disk shows a clear dependence on metal- 
licity, reminding of the oxygen behavior. For argon the 
metallicity effect is less marked, and its evolution remind 
that of intermediate a-elements like silicon or calcium. In 
the absence of stellar observations for these elements, it is 
gratifying that their solar abundances are well reproduced. 

Both odd-Z elements, phosphorus and chlorine, present 
solar or slightly subsolar abundance ratios in the halo- 
thick disk. When the assembling of the thin disk begins, 
their abundances gradually climb towards the solar values 
due to the odd-even effect in massive stars (and in a much 
lesser extent, by contribution of Type la supernovae), until 
the iron from Type la supernovae flattens the evolution of 
phosphorus, although it is not the case for chlorine. 
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Fig. 13. Evolution with metallicity of the iron peak ele- 
ments scandium, titanium and vanadium as a function of 
[Fe/H] . 



4-6. The Iron Peak elements Scandium, Titanium, 
Vanadium, Chromium, Manganese, Cobalt and 
Nickel 

All these elements are basically synthesized through explo- 
sive burning in Type II and/or Type la supernova explo- 
sions, and some of them can also be produced by neutron 
captures during hydrostatic He and C burning. Their stel- 
lar yields are therefore much more uncertain than those 
of the intermediate elements studied above due to our de- 
ficient understanding of the parameters that determine 
supernova explosions (energy and mass cut in Type II su- 
pernovae, progenitors, burning propagation and rates of 
Type la supernovae), and of neutron fluxes in hydrostatic 
He and C combustion. Isotopes with A < 56 are mainly 
produced in explosive O and Si burning, and in Nuclear 
Statistical Equilibrium (NSE), while isotopes with A>56 
are produced in alpha-rich freeze-out NSE, besides contri- 
butions from neutron captures. 

The calculated history of scandium, titanium and 
vanadium is displayed in Fig. 03, and that of chromium, 



manganese, cobalt and nickel appear in Fig. 14 



The evolution for the [Sc/Fe] ratio agrees reasonably 
well with the data for halo stars, characterized by a flat 
and nearly solar value, although between [Fe/H] = -2.5 
and [Fe/H] = -1 we get a gentle decrease and our results 
fall below the observations. At the transition to the thin 



disk the abundance grows a bit, as expected for an odd 
element, but iron from Type la supernovae, which make 
little scandium, overwhelms this increase and lowers the 
ratio [Sc/Fe] during the disk evolution, which is still below 
solar. 

Titanium, dominated by 48 Ti, should behave as an a- 
element. In fact, the observations point to [Ti/Fe]~0.3-0.4 
in the most metal-poor halo stars, and to a clear trend 
towards a declining ratio as metallicity increases. Loosely 
speaking, our calculations reproduce such characteristics, 
but the abundance is systematically lower than observed 
for the complete range of metallicity. At the solar birth all 
the titanium isotopes except 46 Ti are underabundant (see 
Fig. |). TWW |1995| and GP |2000| found the same problems 
with titanium. The conclusion that current yields from 
massive stars alone fail to explain the titanium history is 
unavoidable. 

The isotope 51 V is predominant in the total abundance 
of vanadium. The few available data show a [V/Fe] flat 
and solar in population I and in metal-rich halo stars, and 
a tendency to increase as [Fe/H] decreases from -2 to -3. 
The calculations give a flat [V/Fe] in the halo-thick disk 
with a value of ^-0.3, a moderate and rapid rise at the 
beginning of the thin disk formation, and a flat evolu- 
tion slightly lower than solar during most of the thin disk 
phase; here, again, iron from Type la supernovae compen- 
sates the vanadium produced by massive stars with Z > 
O.1Z . 

The observations on the ratio [Cr/Fe] available until 
mid-nineties showed a constant and solar value down to 
[Fe/H]~-2.5, i.e. the same behavior as iron, indicating that 
both elements were produced in the same scenarios. How- 
ever, recent observations show that the chromium abun- 
dance ratio to iron decreases as the metallicity goes down 
from [Fe/H] =-2.5 to [Fe/H] = -3.5, reaching at the latter 
metallicity values as low a s [Cr/Fe]^-0.7. The chromium 
stellar yields of WW 1995] are quite insensitive to Z. Our 
results, shown in Fig. [l4|, give an almost flat [Cr/Fe] in the 
complete range of metallicities, and thus fail to reproduce 
the data at low metallicity. The flatness of the abundance 
ratio is more marked in the disk because chromium and 
iron are produced in the same proportions by Type la 
supernovae, fitting precisely the solar abundances of the 
chromium isotopes. 

In the case of manganese, with only one stable and 
odd-Z isotope ( 55 Mn), the calculated evolution fits well 
the data for [Fe/H]>-2.5, although in the thin disk there 
is an slight overabundance. The roughly linear increase 
in the disk reveals the odd-even effect for this element. 
At very low metallicities, as for chromium, the calculated 
abundances lie clearly above the data. 

The observational situation for the ratio [Co/Fe] is 
rather puzzling. Cobalt is an odd element and therefore 
one expects its abundance to diminish with decreasing 
metallicity. Instead, the latest data point to a noticeable 
increase of its abundance at the lowest metallicities. We 
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Fig. 14. Evolution of the iron peak elements Cr, Mn, Co 
and Ni as a function of [Fe/H] 



got a rough agreement down to [Fe/H] = —2.5, where 
[Co/Fe] mildly decline when going to lower metallicities, 
although the calculated history appears to lie higher than 
the observations. Below [Fe/H]= —2.5, the abundance ra- 
tio continues its gentle decline, just the opposite trend to 
that observed. 

Nickel abundances in halo and disk stars have been ob- 
tained in a large number of surveys. Most of them show a 
constant and solar value for the ratio [Ni/Fe] at all metal- 
licities, but the most recent observations seem to indicate 
an increasing ratio towards low [Fe/H], although the scat- 
ter is really large. The models give a continuous increase of 
an overabundant [Ni/Fe] for [Fe/H]>-3. We obtain higher 
values than observed due to the combination of a reduc- 



tion of a factor of two in the WW 1995 iron yields from 



massive stars, together with the large amount of 58 Ni in 
the yields of the W7 model of Type la supernovae, the 
latter being responsible for the positive slope in the thin 
disk. 
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Fig. 15. Evolution of copper and zinc as a function of 
[Fe/H] 



4-7. Beyond the iron peak: Copper and Zinc 

These two elements are produced in significant amounts 
at the different burning stages in massive stars. As for the 
iron peak elements, the total ejected amounts of these nu- 
clei are strongly dependent on the characteristics of the 
supernova explosions. Nevertheless, ratios relative to iron 
could be less affected by astrophysical uncertainties. For 



instance, as shown by WW 1995, in the explosion of a so- 



lar metallicity 35 Mq model, variations by a factor 2 in 
the explosion energy give changes in the [Ni/Fe] ratio of 
less than 10%. This argument in favor of the relatively 
robustness of the calculated evolution of the element ra- 
tios relative to iron also applies to copper and zinc. The 
evolution of both elements is presented in Fig. 

Observations for these elements are scarce because 
they do not present enough suitable transitions. Accord- 
ing to the available data, copper behaves like a typical 
odd nucleus, with a marked odd-even effect which gives a 
continuous growth of [Cu/Fe] with the iron contents. Our 
model nicely fits the data. At low metallicities we obtain 
a flat evolution, and in the thin disk the odd-even effect 
is clearly seen. As shown in Fig. pi the solar abundances 
of the two copper isotopes are well adjusted. 

The data for the zinc abundance show that this ele- 
ment evolves as iron, with a constant and solar [Zn/Fe] 
in the complete interval of observed metallicities. Here, 
again, the calculated evolution fits well the data. However, 
this result could be misleading since even if we obtain an 
elemental abundance in almost perfect agreement with the 
solar value, the distribution between the four isotopes of 
zinc (see Fig. ||) does not agree with the solar isotopic ra- 
tios. As mentioned above, we underproduce by ~ a factor 
of 2 the most abundant isotope, 64 Zn, and overproduce 
68 Zn by the same factor. 
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5. Summary and conclusions 

We have analyzed the evolution of all the stable isotopes 
between hydrogen and zinc (lithium, beryllium and boron 
isotopes will be treated in a future paper) by means of 
a detailed chemical evolution code. In the framework of 
standard chemical evolution models only the works by 
TWW |1995| and by GP |200C| have studied such a large num- 
ber of isotopes. In the first paper, the model employed was 
adequate just for the galactic disk evolution, a part f rom 
using the rather old yields from Renzini & Voli ( |l98l[ ) for 
low and intermediate-mass stars. In the second case, even 
if the halo and the disk were correctly treated, only yields 
from massive stars were included 

Our model assumes the formation of the Milky Way in 
two main episodes of exponentially decreasing infall. The 
first episode forms the halo and the thick disk from in- 
fall of primordial extragalactic material with an c-folding 
timcscale of 1 Gyr. Then, the thin disk is assembled by in- 
fall of metal rich extragalactic gas with a global metallicity 
of O.IZq in solar proportions, and with a relatively long 
timescale (7 Gyr at the solar ring). We have also calcu- 
lated models that accrete material of primordial composi- 
tion during the complete evolution. Although the results 
of the two kind of models are similar, we obtain a bet- 
ter agreement with the data for the stellar distribution 
metallicities in the case of the model with enriched in- 
fall. We have employed an up to date IMF (the version of 
Kroupa et al. (1993)), and a form for the SFR (the Dopita 
& Ryder's (1994) law), that relies on observations of the 
properties of star formation in spiral galaxies. 

One important characteristic of our model is that it in- 
corporates metallicity dependent yields for isolated stars 



of 58 Nii in the adopted yields of Type la supernovae 
and of 62 Ni in the yields of exploded massive stars of 



in the whole range of stellar masses. We use the WW 1995 
yields for massive stars and van der Hoek & Groenen 



gen (1997) for low and intermediate mass stars. Following 
TWW1995, in order to obtain a better concordance with 
observations, we have reduced the rather uncertain iron 



yield of WW 1995 by a factor of two. We have not included 



metallicity effects in the yields of Type la supernovae, but 
GP pOOO pointed out that the inclusion of the metallicity 
dependent yields of Iwamoto et al. (199E) hardly affects 
the final results. Finally we have taken into account the 
nucleosynthetic contributions from novae that give impor- 
tant contributions to isotopes such as 13 C, 15 N, and 17 0. 

We have compared in detail the theoretical results with 
a large ensemble of observational data. First of all, in Sec. 
3 we have shown that our results are able to reproduce the 
main characteristics of the solar vicinity, especially the G- 
dwarf distribution of the local disk, that we consider the 
most stringent restriction. We also obtain a really nice 
agreement with the solar composition (both elemental and 
isotopic) at the time of 8.5 Gyr. It is worth to stress that: 

— All the solar elemental abundances are reproduced 
within a factor of two, except for nickel, which is 
slightly oversolar due to the large production factor 



WW 1995 



— The contributions from low and intermediate-mass 
stars suffice to build the abundances of the main CNO 
isotopes. The inclusion of novae increases the solar 
abundances of 13 C (by a factor of ~1.6), and notice- 
ably that of 15 N (by a factor of ~4) and 17 (by a 
factor of ~2). 

— Massive star yields still encounter difficulties to ac- 
count for the underproduction of the neutron-rich iso- 
topes 48 Ca and 50 Ti, as it is also the case with 47 Ti. 

The evolution of the abundance ratios [element/Fe] ob- 
tained in our model has been confronted with a large and 
updated set of observations. The main conclusions that 
result from the comparison are: 

— The evolution of the CNO elements is very well re- 
produced by our model. As it is well known, massive 



stars underproduce carbon and nitrogen (TWW1995 
GP ^OOO ). However, the inclusion of low and intermedi 



ate mass stars is enough to reproduce adequately the 
carbon evolution and its solar abundance. The behav- 
ior of nitrogen is again well reproduced, although some 
primary nitrogen from massive stars could be neces- 
sary during the halo phase. For oxygen we obtain a 
clear deviation of the so-called a-plateau which is also 
present, in a lesser extent, in the rest of the a-elements 
. Our results show a continuous decrease of the ratio 
[O/Fc] from the beginning of the halo-thick disk phase, 



as pointed out first by TWW1995 and then by Chiap- 
pini et al. (1999). The slope is still shallower than in 



the da ta fro m Israelian et al. ( 1988| ) and Boesgaard 
et al. ( 1999 ), b ut th e fit with the more recent data of 
Carreta et al. ( |2000D is excellent. 

As with oxygen, most of the a and a-like elements 
show some deviations of the typical a-plateau, in the 
sense that their abundances ratios to iron continuously 
decrease with increasing metallicity due to the com- 
bined effects of stellar yields depending on mass and 
metallicity and the contribution to the iron contents 
for Type la supernovae, in agreement with observa- 
tions. However, some of these elements are underpro- 



duced at all metallicities by the WW1995 yields (as 
well as by the yields from Limongi et al. (2000), ac- 
cording to GP2000), especially magnesium and tita- 
nium. 

In general, the intermediate odd-Z elements show the 
expected theoretical behaviour as a consequence of the 
marked odd-even effect in the stellar yields. However, 
for several of these elements there exist important dis- 
crepancies with the observational data, in particular 
for sodium and potassium. The same discrepancy ap- 
pears in the case of scandium and vanadium which, 
besides, are underproduced by the model at all metal- 
licities. 
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— Although to zero order the evolution of the iron peak 
elements is consistent with the data, there are no- 
ticeable discrepancies. Chromium and manganese are 
overabundant at low metallicities, which may indicate 
problems with the iron yields for stars of low metal- 
licity, while nickel is overabundant at moderate and 
high mctallicity. This last flaw has to be attributed to 
the nickel yields of model W7 for Type la supernovae. 
In any case, as mentioned above, the calculations of 
massive star explosions are still plagued with severe 
problems. 
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